Stretching Homopolymers 



2.5 



Greg Morrison 1 ' 2 , Changbong Hyeon 3 , N. M. Toan 2 , Bae-Yeun Ha 4 , and D. Thirumalai 

o 

q 1 Department of Physics, University of Maryland at College Park, College 

(N Park, MD 20742 

2 Biophysics Program, IPST, University of Maryland at College Park, 

College Park, MD 20742 
3 Center for Theoretical Biological Physics, University of California at San 

Diego, La Jolla, CA 92093 
CN 4 Department of Physics and Astronomy, University of Waterloo, Waterloo, 

, , Ontario, Canada N2L 3G1 

5 Department of Chemistry and Biochemsitry, University of Maryland at 
College Park, College Park, MD 20742 

S 

Abstract 

Force induced stretching of polymers is important in a variety of contexts. We have 
used theory and simulations to describe the response of homopolymers, with N monomers, 
to force (/) in good and poor solvents. In good solvents and for sufficiently large N 
t-h we show, in accord with scaling predictions, that the mean extension along the / axis 

(Z) ~ / for small /, and (Z) ~ /» (the Pincus regime) for intermediate values of /. The 
theoretical predictions for (Z) as a function of / are in excellent agreement with simulations 
for iV = 100 and 1600. However, even with N = 1600, the expected Pincus regime is not 
observed due to the the breakdown of the assumptions in the blob picture for finite N. 
We predict the Pincus scaling in a good solvent will be observed for N > 10 5 . The force- 
dependent structure factors for a polymer in a poor solvent show that there are a hierarchy 
of structures, depending on the nature of the solvent. For a weakly hydrophobic polymer, 
various structures (ideal conformations, self-avoiding chains, globules, and rods) emerge 
on distinct length scales as / is varied. A strongly hydrophobic polymer remains globular 
as long as / is less than a critical value f c . Above f c , an abrupt first order transition to a 
rod- like structure occurs. Our predictions can be tested using single molecule experiments. 
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I. Introduction 



Single molecule nanomanipulation methods have been used to measure the response of biological 
macromolecules to mechanical force. Such measurements give direct estimates of the elasticity of 
DNA [I], RNA [2], proteins [Hill], and polysaccharides [5]. Although tension-induced stretching 
of RNA [El IE] and proteins [JJ [9] largely depends on the architecture of the folded confor- 
mations [TU], sequence effects [llj make it difficult to unambiguously interpret the measured 
force-extension curves (FECs) in terms of unfolding pathways. In this context, stretching of ho- 
mopolymers by force provides a potentially simpler case for which the FECs can be calculated. 

In a pioneering paper, Pincus [12J considered the strong stretching of homopolymers in a 
good solvent. The strong stretching limit corresponds to a large enough force, /, such that 
N v a < (Z) <C Na, where v = 3/5 is the Flory exponent, N is the number of monomers, a is the 
size of a monomer, and (Z) = (z^—Zo) is the mean tension-induced end-to-end distance (we have 
assumed that / is aligned with the z axis). Pincus showed that the size of the stretched polymer 
should be determined by an interplay between the Flory radius R F = N u a and the tensile 
screening length (or the blob size) [12], £p = ksT/f. When / is small, then x = R F /(,p <^ 1, 
while in the opposite limit, x ^> 1. The scaling assumption is that for arbitrary /, the average 
end-to-end distance can be written as 

(Z) = R F $(R F /£ P ) (1) 

With this assumption, one can anticipate three regimes in the FEC. 

(i) For small /, we expect a linear increase in the extension of the chain with f in the z- 
direction. At low forces, $(x) ~ x, and hence (Z) oc R F x (aj3f). (it) In the strongly stretched 
limit, which arises for intermediate forces, the value of (Z) can be obtained by dividing the chain 
into a sequence of aligned tensile blobs (along the force axis) whose size is £p ~ (fifY 1 [13J. 
The monomers contained within each blob behave as an unperturbed self-avoiding walk. In this 
case £p = (/3/) _1 ~ N£, with iV 6 the number of monomers in a blob. The linear extension of 
the chain is then given by (Z) ~ £ P x N/N b ~ Nifif)*- 1 ~ N(Pf)s. We will refer to this 
intermediate scaling regime as the Pincus regime. It should be stressed that this argument is 
valid only if N 3> (£p/a)^ 3> 1, which may not be satisfied for a stiff polymer, or a flexible 
polymer with small iV (see below). (Hi) For extremely large forces (beyond the strong stretching 
regime), we expect the excluded volume to become irrelevant, as the bonds between monomers 
become fully aligned with the z-axis, and no monomer interacts with any other monomer. The 
FEC in this regime will be model-dependent, with (Z) w Na 2 (3f /3 for an extensible chain, and 
(Z) pa Na for a inextensible chain. We will refer to this behavior as the non-universal regime. 

The Pincus scaling description of the stretching of homopolymers is well known. However, 
as far as we are aware, a microscopic derivation of the FEC anticipated by Pincus has not been 
provided. More importantly, it is unclear how the FEC of polymers with finite iV compares 
with the predictions of the scaling theory. In other words, for finite values of iV (on the order 
of 1000), how pronounced is the Pincus regime? In this paper, we develop a self-consistent, 
variational theory based on the Edwards-Singh method [Hj to determine the average extension 
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of a homopolymer in a good solvent. The theory gives excellent agreement with simulations. 
Surprisingly, neither the theoretical predictions nor simulations display the Pincus regime for 
iV = 100 or N = 1600. We show that this is due to a finite-size effect, and show that the 
Pincus regime emerges only for N > 10 5 . Only when N is large is the concept of the tensile blob 
(with £p ~ aN£) satisfied, where N 3> Nb 3> 1. We also show using theory and simulations 
that the tension induced stretching of homopolymers in a poor solvent exhibits a first order 
transition between an ensemble of collapsed states and rod-like conformations. The nature of 
the transition is dependent on how poor the solvent is, which is measured in terms of the relative 
attraction between the monomers. The theoretical predictions for the poor solvent case are only 
in qualitative agreement with the simulations. Simulations of a polymer in a poor solvent show 
that tension-induced transitions occur via a hierarchy of structures, depending on the solvent 
quality. Force-dependent structure factors show that, for a weakly hydrophobic polymer, the 
transition to the stretched state occurs through a variety of structures, depending on the length 
scale (or the magnitude of the wave vector, q). For a strongly hydrophobic chain, the globule 
to rod transition occurs by a first order transition when / exceeds a critical value. 



II. Polymers Under Tension in a Good Solvent 

Theory 

Extensible Polymer: The Hamiltonian for a self-avoiding polymer chain under tension is taken 
to be 

o rN rN 

PH = — dsr 2 (s)-(3f dsz(s) + A 2 , (2) 
Za Jo Jo 

where / is aligned with the z-axis, (3 = 1/ksT, and 

rN rN 

A 2 = ^ / ds ds' S[r(s)-r(s% (3) 
* Jo Jo 

with v the strength of the self-avoiding interaction, with vq > in a good solvent. To compute 
the force-extension curves (FECs) and compare them to simulations, we use a self-consistent 
variational method, originally proposed by Edwards and Singh [14j. Following the convention 
in single molecule experiments, we use FEC for the extension changes upon application of force. 
However, throughout the paper, we will derive and plot the extension (Z) as a function of /. A 
reference Hamiltonian 

o rN rN 

PHl = 2gJP I dsi2 ( s )-PfJ o rfsi ( s )> ( 4 ) 

is chosen and the parameter / is determined self-consistently. Because we are interested in calcu- 
lating the FECs, the relevant quantity is the dependence of (Z(f)) = (z N — z ) = ds (z(s)) 
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on /, where (• • -)o indicates the Boltzmann- weighted average with respect to [3 Ho- In Appendix 
A, we also consider the square of the transverse fluctuations using the Edwards-Singh method. 
Because it is not possible to compute the exact average (Z) , we calculate the difference between 
(Z)q and (Z) 1 (where (• • -}i is the average with respect to (3Hi), assuming that Ai + A 2 is small, 
with 



Ax 



3 
2^ 



1 



JV 



ds r 2 (s). 



To first order in Ai + A 2 , we obtain 

(Z) - (Z)i = (Z(A 1 + A 2 )H - (Z^A, + A 2 >!. 

A self-consistent equation for I is obtained by insisting that (Z) 
condition 



(5) 



(6) 



(Z(A 1 + A 2 )) 1 = (Z) 1 (A 1 + A 2 ) 1 . 



(Z)i, which leads to the 

(7) 



Throughout this work, we compute averages with respect to Hi, so the subscripts on 
be dropped. The terms involving A 1; A 2 are easily calculated using 



•) will 



(ZAi)-(Z)(Ai) = h(f 



n d(Z) = A 2 (A 2 - ljiVa 2 ^ 



and 



(ZA 2 ) - (Z)(A 2 ) 



N 



ds 



01 



N Q 

ds' 
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-(5[r(s)-r(s')]) 



(9) 



o 



d((3f) 

with (5[r(s) - r(s')}) = (3/27ra 2 / 2 |s - s'|)l exp (-|s - s'\l 2 a 2 (3 2 f 2 /6) (the details of the calcula- 
tions are given in Appendix A). Using Eqs. (]7]), (pj), and ([£]), the self-consistent equation for I 
becomes 



/ 2 -l 



vVN 



du- 



ll 



-Nui 2 >f 2 /e 



6v 



u 

2,„2 , 



+ 



3 



1- 



erf I lip 



5iV 



(10) 



where we have defined the dimensionless excluded volume parameter v = (3/27r)5 t> /a 3 , the 
dimensionless force (p = a(3f , and where erf(x) is the error function. We have also included a 
cutoff, 8, in the integral over u (with u = \s — s'\/N), to account for the finite separation between 
the monomers, which is neglected in the continuum representation of the Hamiltonian in Eq. Q. 
We expect 5 ~ l/N , since the discrete monomers are separated by a distance |r i+1 — ^ la on 
an average in the reference Hamiltonian, (3 Hi. The cutoff is only imposed in theories that have 
a self-energy divergence [TH [151 EE 02] > an d is generally not required if there is no divergence, 
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as is the case here. However, we will see that this cutoff is essential in order to reproduce the 
FECs obtained in simulations. Given a solution I to the SCE, the linear end-to-end distance is 
given by (Z) = Na 2 l 2 (3f/3. 

It is not difficult to show that, as / — ► 0, a solution to Eq. (|7| is I ~ 1$ oc (v 2 N)to, giving the 
expected linear regime, (Z) ~ N^v^a x (a/3f). We immediately see that this gives the correct 
scaling with N and v for low forces, with (Z) ~ (R 2 )/ =0 x (/3f)/3. We also note that, if we set 
5 = 0, we exactly recover (in our notation) the original, tension-free self-consistent equation for 
a self-avoiding chain, I 2 — 1 = y^QN/ir 3 voa 3 /l 3 developed by Edwards and Singh [Hj. 

For intermediate /, we can obtain the correct Pincus scaling for large N. If we assume 
/ rs £q, we find Nl^ip 2 ^> 1 when <p »i afifx ~ N~%v~^, defining the transition force fx into 
the strongly stretched Pincus regime. For / > fa, we can neglect terms on the order iV -1 and 
exp(— Nl 2 cp 2 /6) for large N, and set eri(l(p^N/6) « 1. This gives the approximate SCE 



2 f V 67T 

r - 1 ~ — r; 



- erf (^/¥) 



O^" 1 ). (11) 



With 5 ~ Z/iV, we see that we can neglect the error function in this regime as well if a(3fx ~ 
l 2 tiTo iV _ 2i5 <c 1. If N is sufficiently large to satisfy this requirement, the SCE becomes I 2 — 1 ~ 
v\[&k '/fitp + 0(vN~^). We thus find the approximate solution in the Pincus regime I ~ lp oc 

1 12 

(v/(p)B. For large iV and intermediate forces, we find (Z) oc Nv^fa, as is expected [12]. Note 
that neglecting terms of order v may be valid only for extremely large iV (on the order of 
N ~ 10 5 ). Thus, the onset of the non-linear scaling regime depends on both v and N, as was 
anticipated by Pincus. 



For sufficiently large (p, we can neglect terms of order yj -1 in Eq. (10), to find an extended 
or rod-like solution I w 1% — 1. This root gives (Z) w Na 2 f3f/3, identical to the non-interacting 
average for an extensible chain. This is not surprising; as the tension becomes large, the excluded 
volume interaction is not relevant. We also note that, in this regime, the chain will become 
greatly overextended. As was shown by Pincus, the extension beyond the non-linear regime is 
non-universal and depends on the precise model used for the homopolymer [T2] . 

Inextensible Polymer: Because the extensible polymer can overstretch for large forces, which 
may not occur for real polymers that are linked by covalent bonds with high spring constants, we 
develop a theory for an approximately inextensible model. We were also motivated to consider 
the inextensible model because the Monte Carlo simulations for N = 1600 (see Appendix B) 
were performed for a model in which the distance between successive beads is precisely a. We 
begin with the discrete, non-interacting, spring-like Hamiltonian 

n n 

The average end-to-end distance, as well as fluctuations in the x and z directions are easily 
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computed using this Hamiltonian. Denning X 

(Z) 



x N 



Na 

(Z 2 ) - (Zf 



Xo, we find 

\- r f dx x 2 e~ k{x ~ 1)2/2 cosh(ipx) - - 



Na 2 



(X 2 ) 
Na 2 



AT 
1 

I f(Z) 



dx x 3 e k{x 1)2/2 sinh(v9a;) - 



Na 



cp \Na 



<p\Na 



(13) 
(14) 
(15) 



where Af = J °° dx xe k( - x 1)2/2 sinh^x), and we have used (Z 2 ) = N(z 2 ) + N(N — l)(z n )^ 



in Eq. (14). We approximate the Hamiltonian in Eq. (12) with a continuous chain using an 



Inextensible Gaussian Hamiltonian (IGH) [T8l IT9 



Hj[v{s)\ 



N 



ds[ x*(s) + yZ(s) + 



z 2 (s) 



2a 2 . / "V ®l(k,<p) oil{k,ip) 



N 



/3g(k,(p) I dsz(s) 



(16) 



where oti and 0:3 are the effective spring constants in the longitudinal and transverse directions, 
respectively, and g is an effective tension. The spring constants a± and 03, and the effective 
tension g, are functions of k and </?. Using the IGH, we find 



(Z) 

(z 2 ) - (Z) 2 

(X 2 ) 



NaalPg/3, 

Na 2 a 2 3 /3, 

Na 2 al/3. 



(17) 



Equating the averages in Eq. (17) with those in Eqs. ( 13 )-( 15 ) explicitly gives the desired IGH 



in terms of k and /. The full expression for the a^'s and g are quite lengthy for general k and 
/, and we omit them here. Note that, with an insertion of 5(x — 1) into all integrals in Eqs. 
(13 )-( 15 ), or equivalently, in the limit as k — * 00, we recover the Freely Jointed Chain (FJC) 



averages. In the FJC limit, the expressions for the a,'s and g are quite simple, and we find 
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(p COtll(<y9) — 1 



tt'3 

aj3g 



. M 1 - <^ 2 csch 2 (y9) 



These spring constants, ai(k — > 00,/), were derived by Hatfield and Quake using a different 
method [19] . 

We note that this approximate FJC Hamiltonian gives the simple Gaussian behavior for 
tp — >• 0, whereas in the limit of ip — > 00, we can easily show that the distributions give the 
expected form of P{X) = 8(X), P(Y) = 6(Y), and P{Z) = 5(Z - Na), with X = x N - x , 
and similarly for Y. We therefore expect that the IGH to be an excellent approximation for an 
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inextensible chain in the limits of small and large /, with possible deviations from the correct 
distribution for intermediate /. Because of the more complicated form of the Hamiltonian in 



Eq. (16), exact analytic work is difficult in the inextensible case. We can, however, generate a 
self-consistent equation using Eq. ^ to determine the FEC of a self-avoiding inextensible chain 
in a manner similar to the extensible case. Using the reference Hamiltonian 



and defining 



we can, to first order in A[ IGH ^ + A 2 , develop the self-consistent equation (Z(A ( { GH ^ + A 2 )) = 
(Z) (A^ GH ^ + A 2 ), similar to Eq. (H). The form of the inextensible SCE is similar to that of 



(10), with 



/ 2 - 1 = ^ C du I^V^l7V6 ; (21) 
with 7 = a/3g the dimensionless effective tension. It is possible, albeit complicated, to show 



that the solution to Eq. (21), with k ^> 1, will be divided into approximately the same scaling 
regions as we found in the extensible case. The solutions to the inextensible SCE, determined 
using Eq. (21), are similar to the extensible roots from Eq. ([7]), with significant differences 



in the two models occurring only for ip > 1. Again, the expected Pincus scaling of (Z) ~ fs 
emerges only for very large N. Thus, both for the extensible chain and the IGH with excluded 
volume interactions, the linear behavior and the Pincus regime are obtained. The behavior of 
the FEC in the limit of very large force is clearly model dependent, as predicted by Pincus 
[T2] . The theoretical predictions for the IGH with excluded volume interactions are validated 
by explicit comparison to Monte Carlo simulations (see below). 



Simulations 

Extensible Polymer: In order to determine if the theory accurately predicts the effect of excluded 
volume on a self-avoiding polymer under tension, we have performed Langevin simulations with 
iV = 100 at various stretching forces. To calculate the equilibrium FEC of a self-avoiding 
polymer, we performed low friction Langevin dynamics simulations using the Hamiltonian 

o N-l N-2 AT , v 12 

^ = ^E^- r /- fl2 ) + EE £ Ur-^ -Mzn-z!), (22) 

i=l i=l j=i+2 V I 4+1 ' J l / 

with a = 1, e = 100, and N = 100. We set ksT = 1/(3 = 1 in the simulations. The first term 



in Eq. (22) describes the chain connectivity in the extensible form that, in the continuum limit, 
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apf=<p 



Figure 1: (Z) as a function of ip for varying ip. The dots are the simulation results with N = 100. 
The linear regime corresponds to the dimensionless excluded volume parameter v ~ 58.6. The 
best visual fit (solid line) is obtained with 5 = 1.61 /N. Also shown are the fits with 5 = 
(dotted line). The inset compares the theoretical predictions (solid line) and the simulations 
results (dots) for the effective scaling exponent y. 



becomes 3/2a 2 dsr 2 (s). We model the excluded volume interactions between the monomers 
using a r~ 12 repulsion term (the second term in Eq. (22)). Because of the large e value, the 
summation does not include neighboring monomers (i and i + 1) to avoid excessive repulsive 
forces. The last term in Eq. (22) denotes the potential due to tension acting on the ends of the 
polymer. Thus, this model can be viewed as the discrete representation of the Hamiltonian in 
Eq. @). 

The Langevin equations for each monomer are integrated in the low friction limit, which 
has been shown to accelerate the sampling rate of the conformational space of the polymer [20J . 
The equations of motion are 

OH -» 

mr'i = -(ri -— h f i , (23) 

where m is the mass of the monomer, ( is the friction coefficient, —dH/dYi is the conformation 
force arising from Eq. (22), and Tj is a random force that satisfies the fluctuation-dissipation 
theorem, (Ti(t) ■ Tj(t')) = 6(ksT/h 8(t — t')5ij, where the integration time (h) is discretized. 
The natural time is = (ma 2 /e h ) 1 ' 2 . We chose ( = Om^ 1 and h = 0.002r L . To begin the 
simulations, we generate 200 initial random polymer conformations, and thermally equilibrate 
those structures for 5 x 10 6 h with / = 0. Subsequently, a constant force is applied in the 
^-direction to one end of each polymer, with the other end held in a fixed position. The force 
exerted is increased as f) = 10 _3+0 ' 1 - J ksT/a with j = 1,2,..., 39. The integer j is increased 
every 5 x 10 6 h. For each force step, we neglect the first 2 x 10 6 steps to ensure that the chain 
has equilibrated at f), and collect the statistics of polymer conformations every 10 4 integration 
time steps for the remaining time steps. 
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In order to compare our theory to simulations, we need two fitting parameters, Vq and S. 
We determine Vq by fitting the simulated FEC in the linear, low force regime, and obtain 5 
by a global fit of the theoretical predictions of the FEC to the results of the simulation. The 
scaling laws for the extension as a function of force can not be accurately determined by simply 
fitting a linear [2U E2J E3] or log-log (21] plot of the FEC. Such fits implicitly assume that there 
exists a well-defined scaling regime, where (Z) ~ f y with y constant. In order to determine the 
various scaling regimes of the FEC without imposing such an assumption, we will define the 
force-dependent effective scaling exponent y such that 

In Figure 1, we show the best fit of the theory compared to the simulations for the polymer in a 
good solvent (v > 0). With the choice of v ~ 58.6 and 5 = 1.61/N, the theoretical predictions 
agree well with the simulation data. We note that this gives vq ~ 178a 3 , significantly larger 
than the hard-core second virial coefficient of Vq = Aira 3 /3. It is known from the Edwards-Singh 
calculation [H] (with / = 0) that if higher order terms are included in deriving the self-consistent 
equation (Eq. (J7])), they merely renormalize v without altering the scaling behavior. A similar 
behavior is expected when / ^ 0. As a result of the renormalization of v, we find that the 
extracted value of v from simulations is larger than the naive value calculated from the second 
virial coefficient. 

We see that the theoretical predictions depend very strongly on the choice of cutoff, with 
the 6 = theoretical FEC showing very poor agreement with the simulated data for a(3f > 0.1. 
This is somewhat surprising, as a cutoff in the continuum limit approximation generally is used 
only to avoid self-energy divergences in the theory (TU [151 EH EL which are not present here. 
We also note that neither the theory nor the simulation predicts a Pincus-like scaling of y ~ 2/3, 
because the notion of the unperturbed tensile blob is not applicable for iV = 100 (see below). 

In order to asses the conditions under which the Pincus regime can be obtained, we plot 
the theoretical effective scaling exponent y((p) for increasing iV in Fig. 2 (a). While there is 
no clear Pincus regime for N=100, the expected 2/3 scaling emerges for larger N. Variation 
in v (i.e. changing the interaction strength of the excluded volume) only effects the depth of 
the trough (see Fig. 2(a)) in the final transition (data not shown), so adjusting v can not yield 
the expected Pincus scaling for smaller N. Figure 2 also shows that a very large N ~ 10 6 is 
required in order to see the 2/3 scaling over a large force range. For small values of N, the 
inequality N ^> (£p/a)* ^> 1 required to observe the Pincus scaling is not satisfied. The width 
Af over which the strong stretching is observed can be computed using the self-consistent the- 
ory. If we define the Pincus regime such that dy/dip < e (with y defined in Eq. (24)) for some 
tolerance e, we can numerically determine the dependence of the width of the Pincus regime 
with respect to N. The width of the Pincus regime, Af, is shown in Fig. 2 (b), along with a 
fit Af w 0.018 — 1600A^ _1 . In the inset, we show the transition force into the Pincus regime, 

3 

fx, along with the expected scaling of N~$. We can extrapolate that the minimum number of 
monomers, N min , for a self-avoiding polymer to show that the Pincus regime emerges only when 
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Figure 2: (a) The effective scaling exponent y for iV = 10 2 (— ), N = 10 4 (■••), N = 10 6 ( ), 

N = 10 s ( ), and N = 10 10 ( ), all with v = 58.6 and 5 = 1.61 /N obtained theoretically. 

The inset shows the log-log plot of the extension vs. force, for the same parameters, (b) The 
width of the Pincus regime Af as a function of iV for e = 0.05. The inset shows the initial 

3 

Pincus transition force as a function of N. Also shown is the predicted N~$ scaling. 



N m in ~ 9 x 10 4 for e = 0.05. Larger values of A" are required for the Pincus scaling to continue 
over an observable interval of /. This finite size effect is remarkable, because when / = the 
exponent [y « 0.6) can be accurately obtained with A" < 100 [25J. Because N min is too large for 
accurate simulations, it is not possible to explicitly demonstrate the nonlinear scaling in silico. 
In principle, single molecule AFM or optical tweezer experiments can be used to confirm the 
predictions. 



Inextensible Polymer: In order to test our inextensible theory, we determine the best fit to a 
Monte Carlo simulation of a thick chain [2U 122] with A" = 1600. The thick chain is an inexten- 
sible, hard-core excluded volume model, with a configuration rejected if a triplet of monomers 
lie within a circle of radius a (see Appendix B for details). Our variational Hamiltonian in Eq. 
(19) is generated using the spring constant k = 10 4 in Eq. (16). In Fig. 3, we compare the 
FEC and effective scaling exponent (Eq. (p4|)) for the simulations and the inextensible theory, 
in Eq. (21). The FEC obtained using Monte Carlo simulations is in very good agreement with 
the theoretical predictions (Fig. 3(a)). We find v ~ 15.7 gives a good fit for the simulation data 
for low forces, and again 5 = 1.61 /N gives a good global fit to the simulated data. 

In Fig 3(b), we see that there is a deflection in y m 2/3 at (p ~ 0.1, corresponding to a 
Pincus-like regime observed in the simulations, and predicted by the theory based on Eq. (21). 
Such a deflection near y = 2/3 is predicted by the theory for both the extensible and inextensible 
Hamiltonians, and can be clearly seen in Fig. 2(a) for A" = 10 4 . This deflection shows that 
the Pincus regime is beginning to emerge, but the width of the regime Af is vanishingly small. 
We also see the expected return to the non-interacting FJC behavior for large /. The fit is, 
however, quite poor for af3f ^1 — 4, where the effective scaling exponent differs greatly from the 
simulation data. Figure 3(a) shows that the poor fit for intermediate / originates with a slight 
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Figure 3: (a) Force extension curve for an inextensible chain with iV = 1600. Shown are the 
simulation data (dots), along with the best fit for the IGH (solid line) and extensible Hamil- 
tonian (dashed line), with v = 15.7 and 5 = 1.61 /N. (b): Effective scaling exponent y for the 
inextensible FJC. The solid line shows the theoretical exponents for the IGH and the dashed 
lines correspond to the extensible Hamiltonian. Also shown is the non-interacting FJC exponent 
(dotted line). 



overestimation in (Z) vs. f near a/3f ~ 1, followed by an underestimation in (Z) near a/3f ~ 3. 
This over- and underestimation produces a FEC that is not monotonically increasing with /, 
which is a completely non-physical result. The small differences between the theoretical and 
simulated FECs are greatly exaggerated by the effective scaling exponent in the intermediate 
force range. 

The reason for the discrepancy between theory and simulation for intermediate forces is that, 
in the approximate representation for the (nearly) inextensible chain, extensions from |Ar n | = a 
are allowed (see Appendix B). As a result, the chain can stretch somewhat, with mean monomer 
spacing exceeding a. For this reason, less force is required to extend the chain at intermediate 
forces, producing an overestimate of the FEC. The minor disagreement between the theory and 
simulations in the FEC is amplified when the effective exponent y = <91og( (Z) )/dlog(<p) is com- 
puted (Fig 3(b)). We see, however, that both the extensible and inextensible polymer models 
in a good solvent accurately predict the Pincus-like regime observed for a(3f ~ 10~ 2 — 10 _1 . At 
high forces, the response to the force depends on the precise model used to account for chain 
connectivity. As a result, the predictions for the extensible and inextensible polymer models are 
vastly different when a/3f > 1. 



Reexamination of the Blob Concept for Finite N 

In order to better understand the unexpected scaling behavior of the FECs for finite N, a more 
detailed study of the physical processes of extension are required. There are three mechanisms 
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by which the average extension of an extensible chain can increase as a function of force. The 
first is orientation of the polymer along the force axis. We expect that, for small /, the force 
will cause alignment with the z axis, with little perturbation of the chain conformation. In 
the second mechanism, the extension of the polymer is determined by an interplay between 
£p (a length scale below which / is not relevant), and N (which effectively determines the 
number of aligned blobs along the force direction). We expect this mechanism will occur for 
intermediate forces, and for sufficiently large N, cause the emergence of the Pincus regime. As 
these blobs are stretched, (Z) will increase without significantly affecting the alignment along 
the / axis. For large forces we expect overextension to dominate, when the chain is fully aligned 
and the monomers on a length scale £p are stretched. In order to see these physical mechanisms 
of the extension in the simulations of finite, extensible polymers with N = 100, we compute 
the effective force-induced alignment exponent u, given by (Z/|R|) ~ f u , and the effective 
overextension exponent p, given by (L/N) ~ / M . If the polymer is perfectly aligned along the 
z-axis, we expect that the exponent uo — > 0. The variations of the effective exponents u and p 
for N = 100 as / changes are shown in Fig. 4(a). We see that the polymer aligns with the z 
axis at relatively small forces, with full alignment (uj — > 0) occurring for tp 0.1. Overextension 
does not begin until p ~ 3 (in the non-universal regime, see Fig 4(a)), giving a wide range of 
forces in which stretching of the monomers inside of the blobs contributes to the behavior of 
(Z). Representative snapshots of the chain configuration in the three regimes are shown in Fig. 
4(b). 

The absence of a clear signature of the Pincus regime, even for N = 1600, is intimately 
related to the breakdown of the inequality iV ^> (£p/a)" ^> 1. For large enough N, when the 
nonlinear regime in the FEC is observed (Fig 2(a)), the size of the blob £p ~ /cpT/ f is expected 
to scale as £p ~ aN£, where iV& (presumed to be much greater than unity) is the number of 
monomers inside of the blob. The monomer density, p m , inside the blob will scale as 

Pm ~ N b /e P ~ ej U ' 3 ~ (f/k B T) 3 -^. (25) 

In good solvents, v = 3/5, and hence p m is given by p m ~ f^ 3 = f m . If the effective value for m 
with finite N exceeds m = 4/3, as could be the case when the force locally stretches the chain 
segments inside £p, we will find (Z) ~ f x with x ^ 2/3 in the intermediate force regime. 

In order to provide insights into the effective blob response to / for the self-avoiding extensible 
polymer of the finite size (N=100), we have calculated the dependence of the monomer density 
inside the blob on /. To obtain the scaling behavior between monomer density and the force 
from the simulations, we perform the following steps: 

1. Make a sphere of radius b = £p/2, with £p(= fcpT//), centered on the i th monomer and 
count the number of monomers (iVj,) within the sphere whose volume is b 3 . The density 
of monomers within the sphere center at i th monomer is p m {i) ~ N^/b 3 . 

2. Move to the (i + l) th monomer, and compute the density again. 

3. When i — N, the average density is computed using (p m (f)) = 1/N J2f =1 p m (i)- 
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Figure 4: (a) Contributions to the effective scaling exponent y. Shown are y (dots), as well 
as the alignment exponent u (solid line) and overstretching exponent p, (dashed line), (b) The 
blob density as a function of force. The inset shows a log-log plot of the monomer density inside 
the blob, showing three distinct scaling regimes. Scaling relation p m vs / is obtained by fitting 
the data above red line in the inset, (c) The ensembles of structures at a(3f = (red), 0.5 
(green), and 8.0 (blue) are given to demonstrate the three step mechanism of the extensible 
chain stretching, i.e., (i) alignment, (ii) disruption of tensile blob, and (iii) overstretching. 



4. Repeat this procedure for the ensemble of structures obtained at each force. 

Although this method of computing the monomer density from the polymer structures is very 
crude, the scaling exponent between p m and / should not be affected by the details of the 
calculation. The results are shown in Fig. 4(c). We find that p m ~ J 1 ' 6 in the intermediate 



force regime (data above the red base line in the inset of Fig. 4(c)) . From Eq. (25), a den- 
sity scaling of / 16 implies £p ~ A°' 71 ^ N® m , which indicates that there is no force range in 
which ideal blobs can be observed for small N. In other words, the separation in length scale 
iV ^> Nb ^> 1 is not satisfied. The observed scaling exponent for ]V& is greater than that for a 
simple self-avoiding walk, which suggests that the monomers inside of the blob do not behave as 
unperturbed SAW's. Thus, the fundamental premise used in the blob argument used to derive 
the Pincus regime breaks down for small N. The tensile force is felt by the monomers within the 
blobs, which swell due to the stretching of monomers inside £p. The density of monomers inside 
the blob scales differently than the expected for large values of N, and provides the microscopic 
reason why, in the finite-sized self-avoiding chain, (Z) ~ f x with x < 2/3. As N increases, 

4 

the intermediate force regime can be large enough so that p m ~ fs, which is needed to see the 

2 

Pincus scaling (Z) ~ /s. 
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III. Homopolymer in a Poor Solvent 



Theory 

In a poor solvent, the second virial coefficient (vq) becomes negative. The strength of the attrac- 
tive interactions between the monomers exceed that between the monomers and the solvent. As 
a result, the polymer adopts collapsed, globular conformations at temperatures below the Flory 
6 temperature. In poor solvents, the Edwards model is modified to include an effective three- 
body interaction, to ensure that the averages of physical observables converge. The extensible 
Hamiltonian in a poor solvent is (3Hp = (5Hq + A3, where Hq is defined in Eq. (|2]) and 

A 3 = — [ ds f N ds' f N ds" SMs) - r(s')l S[r(s') - r(s")}. (26) 
6 Jo Jo Jo 

The self-consistent equation for the extension in this case becomes (Z(A 1 +A 2 +A 3 )) = (Z)(A 1 + 
A2 + A3), similar to Eq. ([7]). We have already determined the Ai and A2 terms, and need only 
compute (ZA 3 ) — (Z)(A 3 ) = ad/d(p(A 3 ). The SCE for an extensible polymer in a poor solvent 
can be written as, 



/ 3 .1* \fu 

•1 



+» f iui I '"^ (1 -"': / r ) '" 1+ ^ ) e-"V(».^)/«, (27) 



76 / 1 / z 3/2 3/2 

L -IS J5 U{ U 2 



where we have defined w = (3/27r) 3 w^/a 6 . The inextensible self-consistent equation is similar, 
and has a similar root structure, and we will therefore omit such a calculation here. Again, 
we have included a cutoff in the integrals, as was done for the two-body case. However, the 



three-body integral in Eq. (27) is clearly divergent for 6 = 0, unlike the two-body term. This 
divergence must be removed for the self-consistent equation to converge in the limit of — > 00, 
by renormalizing w. For / = 0, we can evaluate the three-body integral exactly, and find that, 
with 5 ~ l/N and as — > 00, it diverges as ~ 16/3y/5. The — > 00 divergence is therefore 
removed if we renormalize w = w/y/N. It is not difficult to show that the self-consistent 



equation has a solution I w l g = (4w/\v\)$ N^^ for / = and large N, giving the expected 
scaling (R 2 ) ~ A^3 for a homopolymer in a poor solvent. However, the final term of Eq. (27) 
can not be evaluated exactly for non-zero (p, so we must resort to numerical work in order to 
determine the roots for larger forces. 



We find that Eq. (27) has three unique roots beyond a critical force f c , which correspond to 
collapsed (l c ), extended (Ie), and saddle point (Z&) structures. Numerically, we find < l c < l g , 
and Ie ~ 1 for / > f c . Our interpretation of l c as corresponding to a collapsed state is 
only qualitative, because an extensible homopolymer (used as the reference Hamiltonian in the 
calculations) in a poor solvent does not have a unique 'collapsed' state. With the interpretation 
that l c and Ie are the roots signifying the two local minima of the free energy for the collapsed 
and extended states, we can interpret the saddle point solution as a local maximum in the 
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(b) (c) 




Figure 5: (a): The extension (Z)/L = (a/3f)l 2 /3 for the three roots of the self-consistent 
equation in a poor solvent for v = — 5 and w = 1, with 5 = 1.61/N: l g (solid line), h (dashed 
line), and 1$ (dotted line). The three values for I are shown in the inset, (b) For / < f c , 
the polymer is globular, (c) In the force range f c < f < the chain conformations are 
a combination of globular and extended states, (d) at f — the globular configuration is 
marginally stable, (e) For / > the chain is in the fully extended state. 



free energy, i.e. the barrier (or saddle point) between the two states. Again, this interpretation 
is qualitative only, because there is no well defined 'barrier' between the collapsed and extended 
states. In Fig. 5(a) we show the extension (Z)/L = (pl 2 /3 for the three solutions to the self- 
consistent equation (27) for v = —5 and w — 1 (arbitrarily chosen), and with 5 = 1.61/N. We 
see a(3f c ~ 3.5 is the critical force at which the extended and saddle point solutions emerge. 
The critical force f c depends on the particular values of v and w, and we expect it will be 
an increasing function of \v\/w. In this triple-root regime, the polymer will be in bistable 
equilibrium between the collapsed ensemble and extended state, suggesting the development of 
a pearl-necklace structure for intermediate /. The collapsed and saddle point solutions coalesce 
for a finite / = f& (Fig 5(a) inset). For / > l c and Z& vanish, leaving the extended root Ie 
the only solution to Eq. (27). This shows, as expected, that the inter- monomer interactions 
become irrelevant for sufficiently high force, and (Z) ~ Na x (a(3f)/3 as / —>■ oo. Schematic 
pictures of the free energy as a function of the extension (Z) (Fig. 5(b-e)) for varying force 
illustrate our qualitative interpretation of the solutions to the self-consistent equation (27). 

A similar multi-root structure has been previously predicted for a polymer in a poor solvent 
with electrostatic interactions [T71 EZ] • These references note the emergence of multiple roots 
beyond a critical value of the backbone charge density (in this respect, equivalent to the tension), 
and qualitatively identify the meaning of the multiple roots as we have. However, because the 
Edwards Singh method can not predict the barrier height or the depth of the minima, we can 
not quantitatively predict (Z) for a polymer in a poor solvent. The qualitative picture, namely 
the tension-induced globule to rod transition which should occur when / > f c , is confirmed 
using explicit simulations of force- induced stretching of a homopolymer in a poor solvent. The 



15 



simulations (see below) also provide a microscopic picture of the structural transitions that oc- 
cur as w 3 , in Eq. (27), increases. 



Simulations 

The simulation procedure used to study the stretching of a homopolymer in a poor solvent is 
identical to the one described for the good solvent case, except for the Hamiltonian used. The 
Hamiltonian in a poor solvent is 



o N-l N-2 N [ , v 12 / x 6" 

(28) 

where e = 0.5 and 1.5 are used for different solvent conditions, and where the other parameters 
are the same as in the good solvent case. The nature of the polymer is characterized by the 



second virial coefficient v% = j dr{l — e~^ Vint ^}, where Vi nt (r) is the second term of Eq. (28). 
When e ~ 0.3, v 2 approaches zero, and corresponds to the theta condition (T = Tq). By 
decomposing Vi nt (r) into repulsive (v rep (r)) and attractive (v rep (r)) parts of the potential, one 
can write v 2 ~ / dr{l — e~ /3vrep ^ (1 — (3v att (r))} —v (l — Tq/T) where v = / dr (l — e~ l3Vrep ^) 
|28j . Therefore, 

r e «r(i-g. (29) 

We find t>2 = —1.9a 3 (Tq m 1.7 x T, weakly-hydrophobic condition) for e = 0.5 and v 2 = —15.2a 3 
(Tq fa 6.4 x T, hydrophobic condition) for e — 1.5. These estimates for Tq as a function of £ are 
approximate. For our purposes, approximate estimates are sufficient to illustrate the response 
of weakly hydrophobic and strongly hydrophobic chains to force. 

In Figure 6, we show the average linear extension for weakly hydrophobic (a) and strongly 
hydrophobic (b) polymers as a function of force. The weakly hydrophobic polymer does show a 
transition between two linear scaling regimes, with the low force behavior of (Z) ~ (R 2 )/=o x 
{/3f)/3, and the high force behavior returning to the non-interacting (Z) = Na 2 f3f/3. The 
transition is very smooth, and does not show the expected first order transition due to the weak 
nature of the interactions, as shown in the inset. The strongly hydrophobic chain does show a 
first order transition around aj3f c ~ 1.8, but with broad dispersion. Variations in the critical 
unbinding force is substantial from molecule to molecule, due to the microscopic heterogeneity 
of the globular structures. The observed plateau in Fig 6(b) is most likely due to full alignment 
of the globule along the z-axis, as was the case for the self-avoiding polymer (Fig 4(a)), and 
seen in the theoretical predictions (Fig 5). There is a large range of forces over which the FEC 
does not resemble either the globular or fully extended states, showing the bistable equilibrium 
between the two. 
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Figure 6: (a) Extension as a function of force for a weakly hydrophobic polymer (e in Eq. (28) 
is 0.5). (b) Same as (a), except the chain is strongly hydrophobic (e = 1.5). The insets show 
the effective scaling exponent y (Eq. (24)). The transition to the extended state in (a) appears 
continuous. For the strongly hydrophobic polymer, the globule — > rod like transition is sharp. 
The transition force depends on the energetic details of the globule. The heterogeneity of the 
transition is manifested as the broad variations of transition force. The ensemble of structures 
found at the globule-to-rod transition force (f c = 1.8/csT/a) are shown in (c). 



VI. The Scattering Function Under Force 

The analysis using scattering experiments is useful for investigating the overall polymer config- 
urations, because the scattering intensity as a function of momentum transfer (I(q) = (1/N 2 
J2i<j exp (iq • rjj))) provides structural information on all length scales. In contrast, the FEC 
only provides information about the extension of the chain. By comparing with the well-known 
scaling relations of 7(g) with respect to q for various shapes, one can obtain the structures of 
the polymer over all length scales. For example, for the various structures we expect I(q) ~ q~ x , 
with x = 2 (Gaussian chain), x = 1 (rod), x = 4 (globules), and x = 5/3 (polymer in a good 
solvent) [29]. We calculated the scattering intensity by integrating the distance distribution 
function obtained from the ensemble of structures, 

I(q, f)= f d 3 rP(r, f)e^ r = A* f°° dr r 2 P(r, /) (30) 
J Jo Q r 

with q = |q|. In our simulations performed under varying tension values, we obtained 4nr 2 P(r, f) 
directly from the ensemble of structures by collecting the histograms between the interval of 
(r, r + dr) with dr = 0.2a. 

An inspection of the scattering intensity I(q) of a homopolymer in different solvent condi- 
tions, shown in Fig. 7, along with snapshots of representative structures, succinctly summarizes 
the shapes adopted by the polymer as a result of the tension induced structural transitions, (i) 
In good solvents (Fig 7(a)), the entire chain of N=100 is characterized by the tensile blob in 
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Figure 7: (a) /(g) for a homopolymer in good solvent under varying tension, (b) /(g) for weakly 
hydrophobic homopolymer under tension, (c) /(g) for strongly hydrophobic homopolymer under 
tension. The arrows in (a), (b), and (c) indicate increasing / values. The tension-induced 
structural changes of a homopolymer are illustrated in three solvent conditions (good, near 
theta, and poor solvent conditions). 



the absence of force (or for small force), with /(g) ~ g~ 5 / 3 for g ~ 0.1 — I. As / increases, the 
tensile blobs continuously change to the rod state, which is indicated by /(g) ~ g _1 . (ii) For the 
weakly hydrophobic condition (Fig 7(b)), i.e., slightly above the theta temperature, the chain 
displays a hierarchy of structures on distinct length scales. When / is small, both signatures of 
Gaussian coil (/(g) ~ g~ 2 ) and globule structure (/(g) ~ g~ 4 ) are found on small length scales 
g _1 < 1, while the chain is characterized by the polymer in a good solvent for g _1 > 1. As / 
increases, the globule to rod transition of the self-avoiding chain takes place continuously, (hi) 
For the strongly hydrophobic condition (Fig 7(c)), the whole chain is collapsed to compact glob- 
ule (/(g) ~ g~ 4 . The globular structure is maintained so that all /(g) 's are practically identical 
for g _1 > 1 as long as / < f c . When / becomes greater than f c , a sharp transition occurs, 
reflecting the globule (/(g) ~ g -4 ) to the rod (/(g) ~ g _1 ) transition. The first order nature of 
force-induced stretching has been previously described using scaling arguments [30J . 



V. Conclusions 

We have developed a general theory for describing the response of homopolymers to an external 
force for arbitrary values of N, the number of monomers. By using both an extensible and 
inextensible model for the polymer in a good solvent, we show that the theoretical results are 
in accord with the predictions of the Pincus scaling laws. The mean chain extension depends 

2 

linearly on the force for small /, and scales as (Z) ~ /s for intermediate / and sufficiently large 
N. Simulations of an extensible chain with N = 100 and the thick chain model with N = 1600 
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were performed to validate the theory. The theoretical predictions for the force-extension curves 
are in excellent agreement with the simulation results. Surprisingly, the expected Pincus scaling 
is not observed in simulations, even for iV = 1600. The theory predicts that the width Af 
for observing the Pincus regime for N ~ O(10 3 ) is vanishingly small. Only when N exceeds 
~ 10 5 can the strong stretching limit ((Z) ~ /i) be unambiguously observed. The failure to 
observe the Pincus scaling is linked to the breakdown of the notion that the monomers inside 
the well-defined tensile blobs are unperturbed. For N ~ O(10 3 ), the monomers inside each blob 
feel the effect of force, which essentially violates the required inequality iV ^> (£p/a)» ^> 1. 

Applying tension to a polymer in a poor solvent produces a much richer set of structures, be- 
cause of the presence of an additional attractive monomer-monomer energy scale. In the absence 
of force, a polymer in a poor solvent forms a globule at T < Tq. For this case, the theoretical 
analysis predicts that the globule to stretched (i.e. rod-like conformation) transition should 
occur abruptly via a first-order transition when / exceeds a critical force. While the simulation 
results are in accord with the theoretical predictions, they show several structural transforma- 
tions, depending on the quality of the solvent. The hierarchy of structures are reflected in the 
force-dependent structure factor. For weakly hydrophobic polymers (T m Tq) for small forces, 
the scattering function I(q) shows signatures of Gaussian and globular structures at large length 
scales (small q), whereas over small length scales the polymer behaves as a self-avoiding chain. 
At large forces, the transition to a rod-like conformation occurs. These structural transitions 
occur continuously as / increases for a weakly hydrophobic chain. Strongly hydrophobic chains 
(T < Tq) adopt globular structures for small forces. The conformation remains globular as long 
as / < f c ~ kBTe/Rg. The globular nature of the conformation is reflected in the I(q) ~ q~ A 
scaling. If / > f c , there is an abrupt transition to the rod-like state, which is reflected in the 
I(q) ~ q^ 1 scaling. 

The predictions made here can be, in principle, validated with single molecule AFM or optical 
tweezers experiments. Our simulations show that the forces required to stretch the homopoly- 
mer (N « 100) is on the order of about 30 pN, which are easily accessible in current experiments. 



VI. Appendix A: Self-Consistent equation for / 



In this appendix, we provide the details for the calculations of (Z)i and (<5[r(s) — r(s')])i that 
appear in Eqs. ^ and (J9J). 



JVr(s)Ze-^ d 
( Z )l = -T^/XZZW, = 7^ log 



JVr(s)e-^ d(/3f) 



Vr(s)e 







d((3f) 



log 



-Na 2 l 2 j3f 



(31) 
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(5[r(s) - r(s')])i = J J J J dr{N)dr{s')dr{s)dr{0) G(r(N) - r{s')\N - s')5[r{s) - r(s')] 

x G(r(s') - r(s)\s' - s)G(r(s) - r(0)|0) 
= G(0\s'-s), (32) 

where the propagator G(- • • ) is decomposed into transverse and longitudinal components, G(R\N) 
G±(R,±\N)Gh (Z\N; f), each of which can be exactly obtained for the reference Hamiltonian. We 
find 



G ± (R ± \N) 



JVr x (s)6 (/„" dsr ± (s) - R ± ) e'^P ! ° ds *^ s) 

jVr L (s)e-^P ! " ds * 2{s) 
JVr ± (s) J ^ e ik <^ dsrx(s) - R± ) e -^^ Vdsi ' 2(s) 



d 2 k 

3 



3R 2 , 



27rNaWj eXP { 2NaH\ ] ' (33) 
where the Fourier representation of the 5-function used, and then the path integral is evaluated. 
The propagator Gj_(Rj_|iV) is obtained after performing the integral in Fourier space. The 
longitudinal component is also similarly calculated for the Hamiltonian with linear force term, 



G\\(Z\N;f) 



JVz(s)5 (j N dsz(s) - Z) e -^r^W+wr^) 



2ixNa 2 l 2 



1/2 



exp 



2Na 2 l 2 



Z- 



Na 2 V(3f 



Thus, we find 

(5[r( S )-r( S ')])i = G(0| S '- S ) 



2ti\s' — s\a 2 l 2 



3/2 



exp 



s\a 2 l 2 (3 2 f 2 



6 



(34) 



(35) 



In obtaining Eq. ([7]), we used the extension Z as the observable to determine the optimal 
value of I. Alternatively, one can also obtain the SCE for I using the transverse fluctuation of 
polymer R^ = X 2 + Y 2 where X and Y are the projections of the end-to-end distance vector R, 
i.e., (R^_(Ai + A 2 )) — (Rj_)(Ai + A 2 ) = 0. Computations involving R^_ are significantly simpler 
than those involving the end-to-end distance vector, R 2 , because the propagators in the x and 
y directions are decoupled from the force-dependent propagator in the z direction. Using the 
same methods as before with our original variational Hamiltonian in Eq. Q, we find 

2Na 2 l 2 (l 2 - 1) 



(RiA 1 )-(Ri)(A 1 ) 
(RiA 2 )-(Ri)(A 2 ) 



3 

a 2 l 2 v 



pN pN 

ds ds'G(0\s' - s)(s' - s). 
Jo Jo 



(36) 
(37) 
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Using Eqs. (36) and (37) we obtain the SCE for I, 

v\/N 



l 2 -l 



du- 



ll 



-Nul 2 <fi 2 /6 



U 



(38) 



which is identical to the equation obtained using the linear end-to-end distance (Z) as the gen- 
erating observable in Eq. (10). Thus, the computation of the FEC is not dependent on whether 
Z or is used in determining the self-consistent equation. 



VII. Appendix B: The Thick Chain Model 

In order to verify the theoretical predictions for the polymer described by the (nearly) Inextensi- 
ble Gaussian Hamiltonian (IGH) with excluded volume interactions, we have simulated the FEC 
using the thick chain (TC) model for a self-avoiding polymer. In the TC model, the polymer is 
viewed as a chain with a finite uniform thickness D, and is represented as a succession of beads 
with position vectors r , r^. All of the bond vectors Ar„ = r n+ i — r n (n = 0, .., N) have the 



same modulus a. Therefore, unlike the IGH where (|Ar„|) « a in Eq. (16), the bond length 
restriction |Ar n | = a is strictly enforced in the TC model. The interaction potential of the TC 
under tension is given by 

H TC = J2 V ( R ^)-i-(r N -r ), (39) 

i,j,k 

where the first term enforces the self-avoidance, and the second term represents the external 
force. In particular, 

V(R ijk ) = { °' % jk I „ (40) 
J ' \ oo, R ijk < D, 

where R^k is radius of the circle going through the triplet of beads k). Physically, the first 



term in the Hamiltonian (Eq. (39 )) ensures the self-avoidance of the chain by rejecting both local 
self-intersection (the local radius of curvature must be no smaller than D) and interpenetration of 
any two portions of the chain at some finite arc-length. Intuitively, it allows only configurations 
satisfying the thickness constraints, that the radii of circles going through all the triplets of 
beads (i, j, k) be greater than D. 

In order to characterize the stretching response of a thick chain with D/a = 1 and iV = 
1600 we performed Monte Carlo simulations using the following scheme. Starting from an 
arbitrary initial chain conformation that satisfies the thickness constraints, the exploration of 
the available configuration space was performed by distorting conformations by means of pivot 
and crankshaft moves. The new structures were accepted or rejected according to the standard 



Metropolis criterion (the infinite strength of the three-body penalties of Eq. (39 ) was enforced by 



always rejecting configurations violating the circumradius constraints). The relative elongation 
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of the chain was calculated for increasing values of the applied stretching force. For each run, 
after equilibration, we measure the autocorrelation time and sampled a sufficient number of 
independent conformations to achieve a relative error of at most 10~ 3 in the average chain 
elongation. For moderate to high forces, this typically entailed the collection of 10 4 independent 
structures, whereas a 10-fold increase of sampling was required for small forces due to the 
broad distribution of the end-to-end separation along the force direction. For small forces, 
conformational fluctuations can be even larger than the mean extension, which makes achieving 
converged results for (Z) more difficult. 
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